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We present a numerical classification of the spherically symmetric, static solutions to the 
Einstein- Yang-Mills equations with cosmological constant A. We find three qualitatively 
different classes of configurations, where the solutions in each class are characterized by 
the value of A and the number of nodes, n, of the Yang-Mills amplitude. 

For sufficiently small, positive values of the cosmological constant, A < A crit (n), 
the solutions generalize the Bartnik-McKinnon solitons, which are now surrounded by a 
cosmological horizon and approach the deSitter geometry in the asymptotic region. For 
a discrete set of values A reg (n) > A cr i t (n), the solutions are topologically 3-spheres, the 
ground state (n = 1) being the Einstein Universe. In the intermediate region, that is 
for A crit (n) < A < A reg (n), there exists a discrete family of global solutions with horizon 
and "finite size". 



1 On leave of absence from Tbilisi Mathematical Institute, 380093 Tbilisi, Georgia 



1 Introduction 



The interplay of gravity and non-linear field theoretical matter models leads to a wealth 
of new and surprising phenomena. In particular, there has been an increasing interest 
in both the structure and the stability of black hole solutions "with hair". (See, e.g., 
and [0] for some key references.) Moreover, self-gravitating field theories have also 
become very popular in cosmology in connection with various inflationary scenarios, the 
formation of topological defects in cosmological phase transitions, etc. 

In this paper we present and discuss some new solutions with various global prop- 
erties of the Einstein- Yang-Mills (EYM) system with cosmological constant A. For a 
limited range of the "bifurcation parameter" A we find a class of solutions which can 
be viewed as a continuation of the remarkable discrete family of particle-like solutions 
discovered by Bartnik and McKinnon (BK) for A = KJ. m the vicinity of the origin, 
these solutions resemble the BK solitons. However, the solutions are surrounded by a 
cosmological horizon and approach deSitter spacetime in the asymptotic region. For each 
node number, n, these asymptotically deSitter solutions exist only for sufficiently small 
cosmological constants, < A < A CI i t (n), where we determine A cr i t (n) numerically. 

When A exceeds A crit (?7,) we obtain a different class of solutions, for which the 
2-spheres (i.e., orbits belonging to the assumed SO (3) symmetry) reach their maximal 
size outside the cosmological horizon. The position of the maximal sphere, 5^, moves 
inwards as A increases and approaches the horizon when A tends to some special value 
A*(n), say. Outside S^ ax a true singularity develops. This region resembles the interior 
of a black hole solution, whose singularity is also shielded by a horizon. For obvious 
reasons, we call these solutions bag of gold configurations. 

These bag of gold solutions continue to exist for A cr it(n) < A < A*(n), where the 
extremal sphere now lies inside the horizon. An interesting phenomenon occurs when A 
reaches the upper limit A reg (n), for which the singularity approaches the horizon. For 
A = A rcg (n) an everywhere regular, spatially compact solution exists for all n. In the 
special case where n = 1 this is precisely the Einstein Universe with a constant energy 
density of the Yang-Mills field on S 3 . This particular solution has repeatedly been 
rediscovered in the past j3J. For higher node numbers, the spatial part of the manifold 
is a "squashed" 3-sphere, and the solutions can only be constructed numerically. 

As is the case for the BK family, it would be valuable to have an existence proof for 
the compact solutions, probably along similar lines as presented in 0, 0. We would also 
like to mention Ref. [^] on EYM solutions with cosmological constant, which contains 
some partial results of the present paper. 

A crucial issue is the question of stability of the solutions presented in this pa- 
per. However, it turned out that this is a quite involved and subtle problem, mainly 
for topological reasons. We shall therefore present this part of our investigation in an 
accompanying paper ||. 

This article is organized as follows: In the second and third sections we derive the 
basic equations and present some special solutions which can be given in closed form. 
The fourth and fifth sections are devoted to the asymptotically deSitter solutions and 
their analytic extensions, respectively. The bag of gold configurations are described in 
the sixth section. Finally, in the last section, we discuss the globally regular, compact 
solutions. 
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2 Basic Equations 



We consider an EYM model with cosmological constant A and action 
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i * (W-2A) + ^tr(FA*F) 
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where G is Newton's constant and g denotes the gauge coupling constant. Since we 
restrict ourselves to configurations with spherical symmetry, the spacetime manifold 
(M,g) has (locally) the structure of a warped product, M = M x R S 2 , with metric 
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Here, g and g denote the metrics on M and S 2 , respectively, and R is a function on M. 
Throughout this paper, quantities referring to (M, g) are endowed with a tilde and those 
for (S 2 ,g) with a hat. The Einstein tensor for warped product manifolds becomes 

(3) 
(4) 
(5) 

where 1Z denotes the Ricci scalar of (M, g). (Small and capital Latin letters are used for 
indices on (M,g) and (S 2 ,g), respectively; a, b, c = 0, 1 and A, B, C = 2, 3.) With 
respect to the diagonal parametrization of the metric g, 

g = - e 2a ^dt 2 + e 2b ^dp 2 , (6) 

which we shall often use in this paper, the d'Alembertian of a function R, say, and the 
Ricci scalar on (M, g) are 

z a - b R')' 



UR 



and 



n = -2e 



-(a+b) 



(e a - b a'Y 



{e b - a RY 

_6— a 



(7) 



respectively. 

For SU(2), the spherically symmetric gauge potential has the general form 

A = a t p + w [rfid'd + r v sin fidip] + (w — 1) [f^cfa? — sin ^rf^] , (9) 

where a = a^dt + a\dp, and ao, ai, w and are functions on M. Here f p = nV l /2, 
= 9*fp, f^, = c^fp/ sin^ and = (sin^ cos ip, sin^sinyj, cos^), where r % are the Pauli 
matrices. In the static, purely magnetic case the choice a = w = is compatible with 
the field equations. The gauge potential (61) now reduces to 



A = (w — 1) [f<p dd — t# sin i? d(p ] 
In terms of w, the stress-energy tensor has the components 
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with V(w) = (1 — W f 

With respect to the parametrization @ of the metric g, the static field equations 
assume the form 
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where we have introduced e M = R and where eqs. (|T^), ( ]T5| ) and are the |(00 + 11), 
|(00 — 11) and |(00 + 11 — 22 — 33) components of the Einstein equations. We also note 
that the (dimension-full) coupling constant k = 8nG/g 2 can be absorbed by introducing 
the dimensionless quantities Rj a/k, pj a/k and Ak. (We shall often set k = 2 in this 
paper, that is, we measure length, time and mass in units of [Gg 2 c~ A ], [G 1//2 gc -3 ] and 
[g 2 G -1 ], respectively; see 0.) 

We shall use two gauges in this paper, depending on whether or not R has a 
local maximum. Considering solutions for which R has no critical point, we can use 
Schwarzschild coordinates, that is, we are alowed to choose the gauge 

R(p)=P = r. (18) 

It is then also convenient to introduce the functions N(r) and cr(r), defined by 



iV = (dr\dr) = e , a = yj-g = e a+b . (19) 
In terms of this parametrization, the static equations (|T4]), (ISj) and ( ]TT| ) become 
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where a dash denotes the derivative with respect to r. Here we have already used eq. 
( |20| ) in the second and the third equations, in order to eliminate the metric function a . 
The function m(r) is defined by the relation 
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When considering solutions for which R develops a local extremum, we use the 
gauge a + b = 0, that is, we parametrize the static metric by the two functions R(p) and 
Q(p), where 



Q(p) = e 
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The static field equations (|14T), (|T6|) and (|T7|) then assume the form 
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where now Q' = dQ/dp, etc. Using eq. (p5|), the remaining equation (|T5| ) becomes a first 
integral, 



{QR)'B! = k [Qw 



12 



V(w) 
2R 2 



1 - AR 



(28) 



It is clear that this coordinate system is also suited to discuss solutions for which R has no 
critical points. However, for obvious reasons, we prefer to use the familiar parametriza- 
tion (jl~8|), (|I~9|) in those cases. 



3 Special Solutions 

Before we present a classification of the static configurations, we consider some special 
solutions which can be given in closed form. 

First, for R(p) = p and constant Yang-Mills amplitude, we find from eqs. (|25|)-(|2"8|) 
above: 



R(p)=p, w(p) = 0,±1, Q{p) = l- — + n^-±p\ (29) 

p 2p z 6 



with M being a constant of integration. For w = (V = 1) this solution corresponds 
to the Reissner-Nordstrom-deSitter Universe with unit magnetic charge, whereas we 
obtain the Schwarzschild-deSitter solution for w = ±1 (V = 0). 

Next, we consider solutions for which both R(p) and w(p) are constants. For 
V(w) = one easily finds 

R(p) = -j=, w(p) = ±l, Q(p) = -Ap 2 + Ap + B, (30) 

which corresponds to the H 2 x S 2 Nariai solution ||10|| . (Here A and B are constants of 
integration.) If w = 0, we find for sufficiently small values of the cosmological constant, 
A < (2k)~\ 

R 2 {p) = 1±V 2 V 2 A Mp) = o, Q(p) = j- 2 (^- 2 -l) P 2 + Ap + B. (31) 

In the limit of vanishing A the solution with the lower sign reduces to the magnetic 
Robinson-Bertotti Universe (with R 2 = |/t). 

Finally, there exists a solution for which the components of the stress-energy tensor 
assume constant values without w(p) being a constant. This is possible only for the 
special value A = 3/(2/t). In fact, 

R(p) = v^sin(-£=), w(p) = cos(-£=), Q(p) = l (32) 
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describes the static Einstein Universe. 

The above examples indicate that the qualitative behavior of the static solutions 
to eqs. (|T4])-(|iT|) crucially depends on the value of the cosmological constant. In the 
following, we shall present a classification of these solutions in terms of A and the node 
number of w. 



4 Asymptotically deSitter Solutions 

We start our numerical investigation by considering small values of A. For A = the 
regular, asymptotically flat solutions of the EYM equations were found by Bartnik and 
McKinnon in 1988 U and, since then, have been subject to numerous studies (see, e.g. 
0,0,0] and references therein). Each solution has a typical size, R n , where n is the 
number of nodes of the YM amplitude w. In the region R > Rn the energy density of the 
Yang-Mills field decays rapidly, and the metric approaches the vacuum Schwarzschild 
metric. 

For small values of the cosmological constant, A <C ^/R n , the contribution AR 2 to 
the energy density is negligible. For R < R n , one therefore expects that the solutions 
do not considerably deviate from the BK solutions. In the region r > R n , however, the 
effect of A becomes significant, which suggests that the metric approaches the deSitter 
metric. Hence - for sufficiently small values of the cosmological constant - the solutions 
are expected to resemble the regular BK solitons, which are surrounded by a cosmological 
horizon at R ~ 1/ y/X and approach the deSitter geometry in the asymptotic region. 

The numerical analysis of eqs. fl20|) -(|22D confirms these expectations. In order to 
find numerical solutions, we need the formal power series expansions of the equations 
(|2T| ) and (|22| ) in the vicinity of the origin, r = 0, the cosmological horizon, r = r^, and 
for r — > oo. In the vicinity of the origin, the regular solutions behave as follows (k — 2): 

w = 1-br 2 + 0(r A ), N = 1 - (46 2 + -)r 2 + C(r 4 ) . (33) 

3 

Near the horizon, defined by A^(r^) = 0, we find with x = r — r^: 

w =w h + w' h x + 0{x 2 ), N = N' h x + 0{x 2 ) , (34) 



where 



K=Ul-^>-Arl) <0 and «4 = ^ . (35) 
Finally, in the asymptotic regime, r — > oo, we have 

a 3V }W (woo) 1 . 1 , 

N=l- — -±r 2 + [VKo) - |Aa 2 ]l + 0(1) . (36) 

To o V V 

Here, b, r/j, Wh, W^, M and a are six "shooting" parameters. 

In order to obtain numerical solutions to the static equations one starts the inte- 
gration with the expansions ([33|) and (|34]) and tries to match the functions w, w' and 
at some intermediate point between the origin and the horizon. The three matching 
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Figure 1: Asymptotically deSitter solution with A = 3x 10 4 and n — 3. For this solution 
one finds 6 = 0.6998, r h = 98.99, w h = -0.505, = -0.774, M = m(oo) = 0.994, 
a = -37, tr(0) = 2 x 10~ 3 , <r(r h ) = 0.99999 and iV mi „ = 1.7 x 10^ 3 . 
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conditions then fix the values of the parameters b, r h and Wh appearing in eqs. (|33|) 
and fl34D - Subsequently, one uses the remaining three parameters in eq. ([36|) to match 
the solutions obtained from numerically integrating between the horizon and infinity. 
Finally, the remaining metric function a is obtained from eq. (|20"D, where a behaves like 

a = a(O) + 0(r 2 ), a = a(r h ) + 0(r - r h ), a=l + 0{^), (37) 

in the vicinity of the origin, the horizon and infinity, respectively (see Fig.l). 

The numerical procedure yields the following result: For each fixed value of A <C 1 
we recover a family of solutions which correspond to the A = BK solitons. Each 
solution is characterized by the value of A and the number, n, of nodes of w inside the 
cosmo logical horizon. Outside the horizon w tends to a constant value, Woo, say. Since 
Woo 7^ ±1, the YM field gives rise to the magnetic charge V = [2tr(P ■ P)] 1 / 2 = w 2 ^ — 1, 
where 

P = — I F , (38) 
4vr Is* 

and where the integration is performed over the 2-sphere at spatial infinity. Here we have 
used eq. (0) and F = dA + A A A to obtain F = (w 2 — l)f p dQ + (w - l^dw A A. (It is 
worthwhile recalling that the solutions with A = have vanishing magnetic charge.) The 
metric asymptotically approaches the Reissner-Nordstrom-deSitter metric with effective 
charge [P 2 - 2Aa 2 /3] 1 / 2 ; see eq. ([*§). 



5 Analytic Extensions 

In this section we construct the analytic extension for a generic metric of the above type. 
Our first goal is to write the metric g in conformally flat form, such that the spacetime 
metric becomes 

g = a 2 N (dt 2 - d X 2 ) - r 2 dn 2 . (39) 

In order to do so, we need the following essential properties of the solutions discussed 
above: Both N and a are smooth functions, where a(r) is bounded and everywhere 
positive. The metric function N(r) is subject to the boundary conditions N(0) = 1 
and — >• —c 2 r 2 as r — > oo. Moreover, N(r) changes sign exactly once, namely at the 
horizon, A^(r^) = 0. By virtue of these properties, the new radial coordinate Xi 

/ \ f r dr , / \ r°° df . N 

*( r ) = / > r < r h, and x{r) = Xoo ~ / r > r h (40) 

increases from zero to infinity as r runs from zero to r^, and then decreases from infinity 
to Xoo as r grows from to infinity. The constant Xoo is fixed by considering the 
expansion of the above integrals in the vicinity of the horizon, 

X = -7j-hi \r - r h \ + C + 0(r - r h ) , |r-r^|<l, (41) 

and requiring that the constant C has the same value in both cases. Here we have also 
introduced the quantity rj, 

V = - \vN'\ r=rh > 0, (42) 



7 




which does not vanish for a regular horizon. With respect to x> the metric now assumes 
the desired form fl39|) which, in a neighborhood of the horizon, becomes 

g = ±a{r h ) r] e~ 2r,x [1 + 0(e- 2vx )] (dt 2 - d X 2 ) - r W , (43) 

where the plus and minus signs refer to the regions r < and r > r^, respectively. 
Next, we note that C( r ); defined by 

Q(r) = _ e - 2 ^M ; r<r h , and CM = e~ 2vx{r) , r>r h , (44) 

is a monotonically increasing function of r with ((0) = —1, Ci r h) = and £ — > 
exp(— 2i]Xoo) > as r — ► oo. Hence, the inverse function, r(C), is well-defined and 
the function F(Q, 

HO = -^v 2 (r(())N(r(0), (45) 

is therefore smooth and everywhere positive. As usual, one finally passes from the 
coordinates (t,x) to the new coordinates (U, V), where 

U = e~ vx sinh^x, V = e~ nx cosh 77%, r < , 

U = e~ vx cosh 77%, V = e~ vx sinh rjx, r>rh- (46) 

The analytically extended metric eventually becomes 

g = F{Q{dU 2 -dV 2 ) -r 2 (C)dQ 2 , (47) 
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Figure 3: Spacetime diagram for the asymptotically deSitter solution with A = 3 x 10 4 , 
n = 2. 

where ( = (([/, V) = U 2 - V 2 . The two functions F(Q and r(C) can be determined 
numerically (see Fig. 2). For the deSitter solution one easily finds 

HO = ^y 2 , r(C) = C G [-1, 1]. (48) 

The spacetime diagram in coordinates (U, V) is displayed in Fig.3. The spacetime 
manifold corresponds to the region U 2 — V 2 G [— l,Coo]- The qualitative features of the 
diagram are identical with those of the deSitter solution. 



6 Bag of Gold Solutions 

The asymptotically deSitter solutions described above exist only for sufficiently small 
values of the cosmological constant: For each fixed value of the node parameter n, there 
exists a maximal value A cr i t (n), say, beyond which the numerical analysis breaks down. 

Solutions which belong to larger values of A exhibit a local extremum of R and 
cannot be obtained in Schwarzschild coordinates. We therefore pass to a parametrization 
of the metric for which R(p) is a dynamical function and choose the gauge e 2a = e~ 2b = 
Q(p); see eq. ©. 

Equations (p5|)-([28|) yield the formal power series at the origin (re = 2), 

w = i _ bp i + {p% R = p + 0(p 5 ), Q=l + (lb 2 - |)p 2 + 0(p 4 ), (49) 
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Figure 4: The horizon radius Rh versus the cosmological constant A for the n = 1 EYM 
solutions. 

where b is the only free parameter. The numerical integration shows that Q(p) develops 
a zero at some p = ph(b,A), indicating the presence of a horizon. Requiring that all 
curvature invariants remain finite at the horizon yields 



p^ph 



lim JQw' = 0. (50) 



As a consequence of this condition we obtain a family of solutions between the 
origin and the horizon, which are parametrized by a discrete set of values b n (A), where 
n is the number of nodes. The parameters Wh, Rh and R' h entering the power series at 
the horizon, 

w = w h + w' h x + 0(x 2 ), Q = Q' h x + 0(x 2 ), R = R h + R' h x + 0(x 2 ), (51) 

are therefore fixed, once 6„(A) is known. Here, x = p — ph, and Q' h and w' h are given in 
terms of Wh, Rh and R' h : 



Finally, we use this expansions to extend the solution beyond the horizon. The advantage 
of this procedure is that it essentially uses only one shooting parameter, b; the remaining 
parameters are then iteratively determined. 

The numerical analysis reveals the following picture: For each value of the node 
number, the horizon radius Rh decreases monotonically with increasing values of A, 
where Rh — > oo for A — > and Rh — > for A — > A rcg (n) (see Fig. 4). The limiting value 



10 







-1 



- Q(p)/3 




i i | i i i | i i i | i 

A=A i/ 




R(P) 

5' A=A 3 \ 

o 

3 


w(p)\ 




, , , 1 


> i 1 



2 4 6 8 

ln(l+p) 



Figure 5: Change of the topology of the EYM solutions. The solution with A = Ai = 
0.3304 is asymptotically deSitter, whereas the one with A = A 3 = 0.3306 is of the bag 
of gold type. The value A = A 2 = 0.3305 is very close to A crit . The functions Q(p) and 
w(p) for the three solutions are almost identical. 
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Figure 6: The bag of gold solution with A = 0.4 and n — \. 
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A reg (n), for which the horizon shrinks to zero, decreases with growing node number n, 
where A reg (l) = 3/4 and A reg (oo) = 1/4. 

Depending on the position of the maximum of R, one finds three qualitatively 
different classes of solutions, corresponding to the following subdivision of the interval 
(0,A reg (n)) (see Fig.4). 

(A) < A < A crit (ri): These are the deSitter like solutions discussed earlier. The 
function R(p) has no critical points for finite values of p. In the asymptotic regime, 
R'{p) behaves like 

R\p) = R' + O(-) as p^oo, (53) 
P 

where the constant R'^ decreases with growing values of A and vanishes for A crit (n). 
Thus, for A — > A cr i t (n), R{p) develops an "extremum at infinity". The numerical analysis 
(see Fig. 5) suggests that for A = A crit (n), R(p) asymptotically approaches a constant 
value, R(oc) = l/y / A crit (n), whereas w(p) tends to ±1; such that for p — > oo the solutions 
coincide with the Nariai solution (0). The topology of the solutions therefore changes 
for A = A cr it(n) (for the solution with one node one has A cr i t (l) = 0.3305). 

(B) A crit (n) < A < A*(n): For these values of A the function R(p) develops a maximum 
for a finite value p e outside the horizon, p h < p e < oo. Since R" < (see eq. (|25|)), R(p) 
decreases for p > p e and becomes zero at some finite value p S i ng , say. In fact, the metric 
function Q diverges as p — > p S i ng , indicating that the geometry becomes singular. (For 
the solution with one node one finds A*(l) = 0.334.) 

(C) A*(n) < A < A rog (n): The behavior is similar to case (B). However, now R reaches 
the maximal value inside the horizon, p e < p^. Since A reg (n) is the maximal value 
for which the solutions exhibit a horizon, R vanishes still outside the horizon, that is, 
Psing > Ph- Again, Q is unbounded for p = p sing (see Fig.6). 

We call the solutions which exhibit a horizon and for which R develops a second 
zero outside the horizon bag of gold solutions. 



7 Compact Regular Solutions 

Until now we have restricted our attention to solutions which develop a horizon. A 
new and interesting type of solutions is obtained in the limit A — > A reg (ra), where the 
horizon and the singularity merge, ph — > p S ing- In this limit, that is for A = A reg (n), the 
geometry turns out to be everywhere regular, in particular at both zeros of R. Moreover, 
the points where R assumes its maximal value, p e , lies precisely between these zeros and 
the spatial geometry is symmetric with respect to p e . Since, in this case, the manifold 
has the topology of M x S 3 , the zeros of R and the 2-sphere p = p e will be called the 
north pole, the south pole and the equator, respectively. 

For each node number n, there exists precisely one value of the cosmological con- 
stant, A = A reg (n), for which one obtains compact solutions of the above kind. For n = 1, 
the solution is the static Einstein Universe with A reg (l) = 3/(2k), already presented in 
the second section: 

R(p) = \^ sin(-^=), w(p) = cos(-y=), Q(p) = 1. (54) 

The regular solutions with higher node numbers are obtained in the limit A — > 
A reg (n) from the corresponding bag of gold solutions. An alternative method, which 
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Figure 7: The n = 3 compact solution. 



takes advantage of the reflection symmetry with respect to the equator, is to integrate 
the field equations on the "northern hemisphere", say. In order to do so, one has to 
impose boundary conditions at the pole (i.e., the origin, p = 0) and the equator (p = p e ). 
The solutions are then obtained by matching the numerical integrations from the pole 
and the equator. 

The formal power series at the origin involve one "shooting" parameter, b, and were 



given in eq. (f49|) . In order to obtain the series expansions in the vicinity of the equator, 
we have to distinguish two cases: Depending on whether the gauge field amplitude w(p) 
is antisymmetric or symmetric with respect to p e , the regular compact solutions will be 
called odd (w e = 0) or even (w' e = 0), respectively. 

(i) w e = 0: For the odd configurations one finds with x = p — p e (k = 2) 
R = R e + X -R" e x 2 + 0(x 4 ) , Q = Q e + )fi" e x 2 + 0(x 4 ) , 

w = w' e x + 0(x 3 ), (55) 



where the field equations (p5|)-(p8|) imply that w' e , R" and Q" are given in terms of R e 
and Q e , 

w? = ^(^- 2 +Ai? 2 -i), K = -^'i q: = ^-rD- 

(ii) w' e = 0: For solutions with even Yang-Mills amplitude we have 

R = R e + j { Rf ] x 4 + 0(x 6 ) , Q = Q e + ~Q'l x 2 + 0(x 4 ) , 
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Figure 8: The embedding diagrams for the n — 1, 2, 3 compact solutions. 

w = w e + ^w"x 2 + 0(x 4 ) . (56) 

As before, the only free parameters are R e and Q e . In terms of these, w e is determined 
by eq. (f§), 

and u>g, i?^, 4 ^ and Q' e ' are obtained from eqs. (|2o])-(P7D, 

KO 2 = K), = -^«) 2 , q: = |j(2^(^ e ) - Rl) ■ 

In both cases, the free parameters are the position of the horizon, p e , the cosmo- 
logical constant, the shooting parameter at the pole, b, and two independent shooting 
parameters at the equator (for instance R e and Q e ). The values of these quantities 
are presented in Table. 1 for the first five compact solutions. The shape of the metric 
functions and the Yang-Mills amplitude is given in Fig. 7 for the n = 3 solution. 



Tab. 1. Parameters for compact solutions. 



n 


A 


b 


V2p e /n 


Re/V2 


W e 


Qe 


1 


0.75 


0.25 


1 


1 





1 


2 


0.364244 


0.429599 


4.824 


1.0150 


-0.5320 


16.656 


3 


0.293218 


0.508831 


15.63 


1.0757 





88.390 


4 


0.270328 


0.540489 


39.64 


1.0483 


0.23549 


417.12 


5 


0.260895 


0.554021 


88.43 


1.0485 





1409.7 


oo 


0.25 


0.569032 


oo 


1 




oo 
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The geometry of the compact solutions may be illustrated with the help of embed- 
ding diagrams. Consider the 3-dimensional Euclidean space in cylindrical coordinates 
R, Z and ip. A surface S of revolution in this space is characterized by a mapping 
p — ► (R(p), Z(p)), and the induced metric on S is 

9s = (# 2 + Z' 2 ) dp 2 + R(pf d<p 2 . (57) 

On the other hand, the metric of a spacelike section S' (with t = to and •& = 7r/2) through 
the geometry of the compact solutions is given by 

9 S > = 7^)^ 2 + Ripfdv 2 - (58) 

Hence, the two geometries coincide, provided that we choose the function Z(p) according 
to 

Z(p) = Z(0) + [ P y/l-Q(jj)R*(p)-P, (59) 

JO V JO n) 



jQ{p) 

where p G [0, 2p e ] . 

The embedding diagrams (R(p),Z(p)) for the n = 1, the n = 2 and the n = 
3 compact regular solutions are presented in Fig. 8. For n = 1, we obtain the circle 
(R(p),Z(p)) = T/K(sm(p/y/^,cos(p/y/~R)), reflecting the fact that the manifold in this 
case is precisely M x S 3 . The spatial sections of the solutions with higher values of n 
resemble prolate ellipsoids (or "cigars"). 

It is also instructive to draw the embedding diagrams for the solutions with horizon. 
In this case, the domain of integration in eq. (|59|) is p e [0^], which yields half of the 
diagram. At the horizon one has Q = and therefore dR/dZ = 0. Since R{ph) 7^ 0, the 
horizon corresponds to the "throat" of the geometry, which connects the two identical 
patches of the manifold (see the conformal diagram in Fig. 3). The resulting diagrams 
for several n — 1 solutions are presented in Fig.9. The diagrams show that the throat 
becomes narrower as A tends to the critical value A reg , where the manifold splits into 
two separate pieces. 



8 Concluding Remarks 

The features of the static, spherically symmetric solutions to the EYM equations depend 
critically on the value of the cosmological constant A. For every node number n, there 
exists a globally regular, compact solution with A = A reg (n). For A cr i t (n) < A < A reg (n), 
the configurations have "finite size" and exhibit a horizon. Finally, for sufficiently small 
values of the cosmological constant, A < A crit (?2), the solutions generalize the BK solitons 
surrounded by a cosmological horizon. 

In this paper, we have restricted our attention to solutions with a regular center. 
The extension to configurations with an event horizon is expected to be straightforward. 
In fact, we have no reasons to doubt that one will find a similar classification for these 
black hole solutions. 

No globally regular solutions seem to exist for A > A reg (n). In this case, the metric 
function Q(p) is everywhere positive and diverges as p — ► p S i ng , where p sing is the position 
of the second zero of R{p). Such solutions may therefore be considered as bag of gold 
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Figure 9: The embedding diagrams for the asymptotically deSitter solution with A = 0.3, 
n = 1, and for the n = 1 bag of gold solutions with A = 0.6 and A = 0.745. 

configurations without horizon. When A is small and negative (|A| <C 1), the solutions 
resemble again the BK solitons, however, they approach the anti-deSitter geometry in 
the asymptotic region. 

We have also investigated the stability properties of the solutions presented in this 
paper. The stability analysis for these - asymptotically not flat - solutions is, however, 
rather involved. In particular, the fact that the size R of the 2-spheres develops a local 
maximum gives rise to the following difficulty: Either the pulsation equations assume the 
form of a regular, formally self-adjoint system with unphysical degrees of freedom, or one 
isolates the unphysical modes and obtains a singular pulsation equation. The methods 
by which these problems can be solved are presented in an accompanying paper ||, and 
here we merely mention the result: all of the solutions described above are unstable. 
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